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Abstract 

Let g: M — )• C be a C°°-function with all derivatives bounded and let tr„ denote the 
normalized trace on the n x n matrices. In the paper |EM ] Ercolani and McLaughlin 
established asymptotic expansions of the mean value E{tr„(5((X„))} for a rather general 
class of random matrices X„, including the Gaussian Unitary Ensemble (GUE). Using 
an analytical approach, we provide in the present paper an alternative proof of this 
asymptotic expansion in the GUE case. Specifically we derive for a GUE random 
matrix that 

E{tr„(5(X„))} = ^ r g{x)V4-x^dx + ^ ^ + 0{n-^'^-^), 

where k is an arbitrary positive integer. Considered as mappings of g, we determine 
the coefficients aj{g), j G N, as distributions (in the sense of L. Schwarts). We derive 
a similar asymptotic expansion for the covariance Cov{Tr„[/(X„)], Tr„[(7(X„)]}, where 
/ is a function of the same kind as g, and Tr„ = ntr„. Special focus is drawn to the 
case where g{x) = and f{x) = -jj^ for A,/x in C \ M. In this case the mean and 
covariance considered above correspond to, respectively, the one- and two-dimensional 
Cauchy (or Stieltjes) transform of the GUE(n, ^). 



1 Introduction 

Since the groundbreaking paper |Vo] by Voiculescu, the asymptotics for families of large, 
independent GUE random matrices has become an important tool in the theory of opera- 
tor algebras. In the paper [HT2j it was established that if X^""^ , • • • , Xr"^ are independent 
GUE(?2, -) random matrices (see Definition 12 . II below) . then with probability one we have for 
any polynomial p in r non-commuting variables that 

lim ||p(xi"V..,X("))|| = |b(xi,...,a:,)||, (1.1) 

where {xi, . . . ,Xr} is a free semi-circular family of selfadjoint operators in a C*-probability 
space {A,t) (see [VDNj for definitions), and where || ■ || denotes the operator norm. This 
result leads in particular to the fact that there are non-invertible elements in the extension 
semi-group of the reduced C*-algebra associated to the free group on r generators (see [HT2j ) . 



A key step in the proof of fll.ip was to estabhsh precise estimates of the expectation 
and variance of tTn[g{p{x["\ . . . , Xr"'*))], where tr„ denotes the normahzed trace, g is a C°°- 
function with compact support, and where we assume now that p is a selfadjoint polynomial. 
In fact it was established in |HT2] and |HST] that in this setup we have the estimates: 

E{tr„ [g{p{x[''\ . . . , } = T[g{p{x,, . . . , x.))] + 0{n-'), (1.2) 

V{tr4^?(p(xi"\...,X(")))]} =0(n-2). (1.3) 

Furthermore, if the derivative g' vanishes on the spectrum of the operator p{xi, . . . , Xr), then 
we actually have that 

V{tr4^?(p(X;"\...,X(")))]}=0(n-'^). 

If we assume instead that g is a polynomial, then the left hand sides of (11. 2p and (II. 3p may 
actually be expanded as polynomials in n~^. More precisely it was proved in |Th] that for 
any function w: {1, 2, . . . ,p} — )■ {1, 2, . . . , r} we have thal0 

E{tr. [Xir^^xi")) . . . } = Y: (1-4) 

where T(w) is a certain class of permutations of {1, 2, . . . ,p}, and (7(7) G No for all 7 in T{w) 
(see |Thj or |MN] for details). It was established furthermore in |MN] that for two functions 
w: {1,2, ... ,p} {1, 2, . . . , r} and v. {1, 2, . . . , g} — )■ {1, 2, . . . , r} we have that 

E{tr4xi%x5))...x5ijtr4xi;')xi;V--xS E -"'^^^^ (1-5) 

-y£T{w,v) 

where now T{w, v) is a certain class of permutations of {1, 2, . . . , p + g} and again (7(7) G Nq 
for all 7 in T{w,v) (see |MN] for details). 

In view of (ll.4p and (II. 5p it is natural to ask whether the left hand sides of (II. 2p and (II. Sp 
may in general be expanded as "power series" in when g is, say, a compactly supported 
C°°-function. In the case r = 1, this question was answered affirmatively by Ercolani and 
McLaughlin (see |EM^ Theorem 1.4]) for a more general class of random matrices than the 
GUE. More precisely, Ercolani and McLaughlin established for a single matrix X„ (from the 
considered class of random matrices) and any C°°-function g with at most polynomial growth 
the existence of a sequence («j(5'))jeNo of complex numbers, such that for any n in N and k 
in No, 

E{tr„(,(X„))} = ^^ + 0(n-^-). (1.6) 
i=o 

Their proof is rather involved and is based on Riemann-Hilbert techniques developed by 
Deift, McLaughlin and co-authors. In this paper we provide an alternative proof for (11.61) 
in the case where X„ is a GUE(?7., ^) random matrix. For technical ease, we only establish 
(II. 6p for functions in the class C^(]R) consisting of all C°°-functions (7: M — )■ C, such that 

^When r — 1, formula (|1.4p corresponds to the Harer-Zagier recursion formulas (see |HT1) ). 
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all derivatives g^''\ k G Nq, are bounded on M. However, all (relevant) results of the present 
paper can easily be extended to all C°°-functions with at most polynomial growth. For each 
j in N we show that the coefficient <yj{g) is explicitly given in the form: 
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for a certain linear operator T: C^(]R) — t- C^(]R) (see Theorem 13.51 and Corollary 13.61) . and 
we describe aj explicitly as a distribution (in the sense of L. Schwarts) in terms of Chebychev 
polynomials (cf. Corollary 14. 6p . The proof of (11. 6p is based on the fact, proved by Gotze and 
Tikhomirov in [GTl], that the spectral density hn of a GUE(n, -) random matrix satisfies 
the following third order differential equation: 

^K{x) + (4 - x^)h'^{x) + xK{x) = 0, (x G M). (1.7) 

In the special case where g{x) = for some non-real complex number A, the integral 
f^g{x)hn{x) dx is the Cauchy (or Stieltjes) transform G'„(A) for the measure hn{x)dx, and 
asymptotic expansions like (11. 6p appeared already in the paper |APSj for a rather general 
class of random matrices (including the CUE). In the CUE case, our analytical approach 
leads to the following explicit expansion (see Section H]): 

G„(A) = ,„(A) + + !?a(^ + . . . + + 0(n---^), (1.8) 

where 

3j-l 
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^o(A) = ^-i(A2-4)V^ and VjW = C,A>^' - 4)-^-'/' (j G N). 



r=2j 

The constants Cj^r, 2j < r < 3j — 1, appearing above are positive numbers for which we 
provide recursion formulas (see Proposition 14. 5p . 

As for the "power series expansion" of (11.31) . consider again for each n a single GUE(n, -) 
random matrix. For any functions f,g from C^{M.), we establish in Section O the existence 
of a sequence {f3j{f, g))j^f^g of complex numbers, such that for any k in No and n in N, 

k 

Cov{Tr„[/(X„)],Tr„[<7(X„)]} = ^^^ + 0(n-2^-2), (1.9) 

j=0 ^ 

where Tr„ denotes the un- normalized trace on M„(C), and where the covariance Cov[y, Z] of 
two complex valued square integrable random variables is defined by 

Coy{Y,Z} = E{{Y -K{Y}){Z -E{Z})}. 

The proof of (II. 9p is based on the following formula, essentially due to Pastur and Scherbina 
(see fPSl): 



Cov{Tr,„[/(X„)],Tr„[^(XO]}= / (^^^^^— ^) (^^^^^)p.(x, i/) dx dy, (1.10) 

JiR2 V X — y / \ X — y / 
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where the kernel p„ is given by 

Pn{x,y) = J[M^/lx)fn^li^/^y)-Vn-li^/^x)ifn{^/^y)]' 



with (fn the n'th Hermite function (see formula ( 12 .ip below). The essential step then is to 
establish the formula (see Theorem 15. 4p : 

p„(x,y) = ^[JiMhniy) - ^h'^{x)h'^{y) - ^h:ix)h:{y)], {{x,y) G R'), (1.11) 

where hn{x) = hn{x) — xh'^{x), and /i„ is as before the spectral density of GUE(n, ^). Using 
fll.lOP - fll.lip and Fubini's Theorem, the expansion (11. 9p may be derived from (11.60 . 
In the particular case where 

f{x) = , and g{x) = , {x G M) 

X — X fJL — X 

we obtain in Section [6] the following specific expansion for the two-dimensional Cauchy- 
transform: 

k 

Cov{Tr.[(Al - X„)-i], Tr„[(/xl - XJ-^]} = \ ^^^^ + 0{n-''^'% (1.12) 



where the coefficients rj(A,/i) are given explicitly in terms of the functions rji appearing in 

ro(A,M) 

2(A-/.)2 

.(A - x)-' - (A - yr\fifi - x)-i -ill- y)-^- 



(11. 8p (see Corollary I6.3p . The leading term Jj°^^'^2 may also be identified as the integral 



X — y / \ X — y 



p{x,y) dx dy, 



where p{x,y) dxdy is the weak limit of the measures pn{x,y) dxdy as n — )• oo (cf. (ll.lOp ). 
The limiting density p is explicitly given by 

1 4 — xy 

P(a;, y) = . 2 /. 9 /. l(-2,2)(a:)l(-2,2)(y), {x, y G M), (1.13) 
47r -y/4 _ a;2^4 _ 

and we provide a proof of this fact at the end of Section [51 

In the paper |APS] the authors derive for a rather general class of random matrices an 
expansion for the two-dimensional Cauchy transform in the form: 

k 

Cov{Tr„[(Al - X^)-\ Tr„[(/il - X„)-i]} = J] rf,- „(A, p)n-i + o{n-^), 

j=0 

where the leading coefficient do,n is given explicitly by 

4,„(A,.)^^(-,=i|lJ=-l). (1.14, 
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with a the variance of the relevant hmiting semi-circle distribution. In the GUE set-up 
considered in the present paper, a = 2, and in this case it is easily checked that (io,n is 
identical to the leading coefficient 2{x^'^)2 P-l^P - 

The density p given by f ll.lSp has previously appeared in the paper [C-Dj . There the 
author proves that if X„ is a GUE(n, ^) random matrix, then for any polynomial /, such that 
JI2 /(^)"\/4 — dt = 0, the random variable Tr„(/(Xn)) converges, as n — )• 00, in distribution 
to the Gaussian distribution A^(0, o"^), where the limiting variance cr^ is given by 



The density p has also been identified in the physics literature as the (leading term for 
the) correlation function of the formal level density for the GUE (see [KKP] and references 
therein) . 

In a forthcoming paper (under preparation) we estabhsh results similar to those obtained 
in the present paper for random matrices of Wishart type. 

2 Auxiliary differential equations 

In this section we consider two differential equations, both of which play a crucial role in the 
definition of the operator T introduced in Section |3l The former is a third order differential 
equation for the spectral density of the GUE. We start thus by reviewing the GUE and its 
spectral distribution. 

Consider a random n x n matrix Y = {yij)i<i,j<n defined on some probability space 
{Q,3^,P). The distribution of Y is then the probability measure Py on the set M„(C) of 
n X n-matrices (equipped with Borel-cr-algebra) given by 



for any Borel-subset B of M„(C). 

Throughout the paper we focus exclusively on the Gaussian Unitary Ensemble (GUE), 
which is the class of random matrices defined as follows: 

2.1 Definition. Let n be a positive integer and a positive real number. By GUE(n, cr^) 
we then denote the distribution of a random n x n matrix X = {xij)i<ij<n (defined on some 
probability space) satisfying the following four conditions: 

1. For any i,j in {1, 2, ... , n}, Xij = xji. 

2. The random variables Xij, 1 < i < j < n, are independent. 

3. If I < i < j < n, then Re(xij), Im(xy) are i.i.d. with distribution A^(0, |o"^). 

4. For any i in {1, ... , n}, X liii IS 9j real- valued random variable with distribution A^(0, a^). 




Py{B) = P{Y e B) 
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We recall now the specific form of the spectral distribution of a GUE random matrix. Let 
(fo, ipi,ip2, . . . , be the sequence of Hermite functions, i.e., 

Mt) = (2'A;!v^)"'/'iffc(t)exp(-tV2), (2.1) 
where Hq, Hi, H2, . . . , is the sequence of Hermite polynomials, i.e., 

r d'' 1 

H,{t) = {-l)'^exp{t')[—exp{-t')\. (2.2) 

Recall then (see e.g. |HTlt Corollary 1.6]) that the spectral distribution of a random matrix 
X from GUE(n, -) has density 



n—l 

h^{t) = (2.3) 



fc=i 



w.r.t. Lebesgue measure. More precisely, 

E{tr„((7(X))} = / g{t)hn{t)dt, 



for any Borel function (7 : M — M, for which the integral on the right hand side is well-defined. 
Gotze and Tikhomirov established the following third order differential equation for hn- 

2.2 Proposition ( |GT1] ). For each n in N, /i„ is a solution to the differential equation: 

1 



n 



-^h':{t) + {A-e)h'St)+tK{t) = Q, (te 



Proof. See [UTTl Lemma 2.1]. ■ 

Via the differential equation in Proposition 12.21 and integration by parts, we are lead (see 
the proof of Theorem 13.51 below) to consider the following differential equation: 

(f-A)f'{t) + 3tfit) = g{t) (2.4) 

for suitable given C°°-functions g. The same differential equation was studied by Gotze and 
Tikhomorov in |GT2l lemma 3.1] for a different class of functions g in connection with their 
Stein's Method approach to Wigner's semicircle law. 

2.3 Proposition. For any -function (7: M — t- C, the differential equation 

it" - A) fit) + 3t/(t) = git), (t G M), (2.5) 

has unique C°° -solutions on (—00, 2) and on (—2, 00). Furthermore, there is a -solution 
to fl2.5p on all of M, if and only if g satisfies 



j git)V4:-t^ dt = 0. 



(2.6) 
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Proof. We note first that by splitting / and g in their real and imaginary parts, we may 
assume that they are both real-valued functions. 

Uniqueness: By linearity it suffices to prove uniqueness in the case g = 0, i.e., that is 
the only solution to the homogeneous equation: 

{f-A)f'it) + 3tfit)=0 (2.7) 

on (— oo, 2) and on (—2, oo). By standard methods we can solve (12. 7p on each of the intervals 
(— oo, —2), (—2, 2) and (2, oo). We find thus that any solution to (12. 7p must satisfy that 

■ci(t2-4)-3/2, ift<-2, 
f{t) = {c2{A-t'yy\ iftG(-2,2), 

C3(t2_4)-3/2^ ift>2, 

for suitable constants ci,C2,C3 in M. Since a solution to (12.71) on (— oo, 2) is continuous at 
t = —2, it follows that for such a solution we must have ci = C2 = 0. Similarly, is the only 
solution to (12. 7p on (—2, oo). 

Existence: The existence part is divided into three steps. 

Step I. We start by finding the solution to (12. 5p on (— 2, oo). By standard methods, it 
follows that the solution to (12. 5p on (2, oo) is given by 



/,(t) = (t2 _ 4)-3/2 j (^2 _ 4)1/2^(5) + c(t2 _ 4)-3/2^ ^ (3, Oo), C G C), 

whereas the solution to (12.51) on (—2, 2) is given by 

m = (4 - t2)-3/2 ^'(4 _ s^y/'gis) ds + c(4 - t')-'/\ (t e (-2, 2), c G C). 
Now consider the function /: (—2, oo) — R given by 

r (4 _ t2)-3/2 J^2^4 _ ,2)1/2^^^^ if ^ ^ (_2, 2), 

f{t) = llg{2), ift = 2, (2.8) 

[ (t2 _ 4)-3/2 jt^^2 _ 4)i/2^(s) ds if t G (2, oo). 

We claim that / is a C°°-function on (—2, oo). Once this has been verified, / is automatically 
a solution to (12. 5 p on all of (—2, oo) (by continuity at t = 2). To see that / is a C°°-function 
on (— 2,oo), it suffices to show that / is C°° on (0,oo), and for this we use the following 
change of variables: 

y = t^-4, i.e., t=v/^T4, (t > 2, y > 0). 
For y in (0, oo), we have 

f{V^)=y-'^' l^\s'-AY"g{s)ds = y-'" fu"'- ^^y^ du. 

Ji Jo 2y/u + 4: 
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Using then the change of variables 

u = vy, V e [0, 1], 

we find that 



for any y in (0, oo). Now, consider the function 

i{y)= r.V2.^(v^)d,, 

which is well-defined on (—4, oo). By the usual theorem on differentiation under the integral 
sign (see e.g. [Scl Theorem 11.5]), it follows that / is a C°°-function on (—4 + e,K), for any 
positive numbers e and K such that < e < K. Hence / is a C°°-function on all of (—4, oo). 
Note also that 

m = \g{2) = m 

Furthermore, by performing change of variables as above in the reversed order, we find for 
any y in (—4, 0) that 

Ky) = i-y)-''^' f (4-.^)V2^(.)d. = /(V4T^). 



Hence, we have established that f^^J-i + y) = l{y) for any y in (— 4, oo). Since I is a C°°- 
function on (— 4, oo), and since f{t) = l{t^ — 4) for all t in (0,oo), it follows that / G 
C~((0,oo)), as desired. 

Step II. Next, we find the solution to (12. 5p on (— oo,2). For this, consider the differential 
equation: 

{f _ 4)^'(t) + 3t^(t) = (t e (_2, oo)). (2.9) 

From what we established in Step I, it follows that (12. 9 p has a unique solution ip in 
C~((-2,oo)). Then put 

/i(t) = -^(-t), (tG(-oo,2)), (2.10) 

and note that /i G C°^((— oo, 2)), which satisfies (12. 5p on (— oo,2). 

Step III. It remains to verify that the solutions / and /i, found in Steps I and II above, 
coincide on (—2,2), if and only if equation (12. 6p holds. With as in Step II, note that ip is 
given by the right hand side of (12. 8p . if g{s) is replaced by g{—s). Thus, for any t in (—2, 2), 
we have that 

fit)- flit) = fit) +^{-t) 

= (4 - t2)-3/2 ^'(4 _ sY'gis) ds + (4 - t2)-3/2 _ s2)l/2^(_,) 

= (4 - t2)-3/2 _ sy/'gis) ds + (4 - t2)-3/2 r (4 _ sY^^g{s) ds 
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(4-t2)-3/2 y'(4_52)l/2^(3) ds. 



from which the assertion follow readily. ■ 

2.4 Proposition. For any function C^-function (7: M — t- C, there is a unique C°°-function 
f:R^C, such that 



g{t) = ^ £^ g{s)Vi^ds + {t' - 4)/'(t) + 3tf{t), (t G M). (2.11) 



f2 

Ifge C^{R), then f G C^{R) too. 

Proof. Let g he a function from C°^(]R), and consider the function 



1 

3'= = J g{s)V4: - s2 ds. 



Since fl'c(s)"\/4 — ds = 0, it follows immediately from Proposition 12.31 that there is a 
unique C°°-solution / to (12. lip . Moreover (cf. the proof of Proposition 12. 3p . / satisfies that 



m 



■ (f _ 4)-3/2 j^t(,2 _ 4)V2^^(,) d,, if t G (2, 00), 

-{f - 4)-3/2 (s2 - 4)1/2 (7,(-s) if ^ ^ _2). 



Assume now that g (and hence (7c) is in C^(M), and choose a number R in (0, 00), such that 
|5'c(^)| < -R for all t in R. Then, for any t in (2, 00), we find that 

|/(t)| < (t^ -4)-=^/2^ [\s^Y/^ds = ii?(t2 -4)-V2, 



and thus / is bounded on, say, (3, 00). It follows similarly that / is bounded on, say, 
(—00, —3). Hence, since / is continuous, / is bounded on all of M. 
Taking first derivatives in (12.111) . we note next that 

it' - A)nt) + 5t/(t) + 3/(t) = g'{t), (t G M), 

and by induction we find that in general 

(t^ - 4)/('=+i)(t) + {2k + 3)t/('=)(^) + Kk + 2)/('=-^H^) = a^'^Kt). {ken, te M). 

Thus, for t in R\ {-2,2}, 
and 

Since / and gc are bounded, it follows from (I2.12p that /' is bounded on, say, M \ [—3, 3] 
and hence on all of M. Continuing by induction, it follows similarly from (I2.13P that f^'^^ is 
bounded for all in N. ■ 
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3 Asymptotic expansion for expectations of traces 

In this section we establish the asymptotic expansion i \1.6\i . We start by introducing the 
vector space of infinitely often bounded differentiable functions. 

3.1 Definition. By C^(M) we denote the vector space of C°°-functions / : M — )■ C, satisfying 
that 



< oo. 



VA; G No : sup 

Moreover, we introduce a sequence || • of norms on C^(]R) as follows: 

\\g\\^ = snp\gix)l (geC^iR)), 

and for any A; in No : 

||(7||(fc)= max Wg^^^lU igeC^{R)), 

j=0,...,k 

where g^^^ denotes the j'th derivative of g. Equipped with the sequence (|| ■ ||(A;))fceN of norms, 
C^(]R) becomes a Frechet space (see e.g. Theorem 1.37 and Remark 1.38(c) in (Ruj ) . 

The following lemma is well-known, but for the reader's convenience we include a proof. 

3.2 Lemma. Consider C^(M) as a Fechet space as described in Definition \3.1[ Then a 
linear mapping L : C^{M.) — > C^(M) is continuous, if and only if the following condition is 
satisfied: 

VA; G N 3mfc G N 30^ > V^? G C,°^(M): \\Lg\\(^k) < CMkm,)- (3.1) 

Proof. A sequence ((?„) from C^(M) converges to a function g in C^(R) in the described 
Frechet topology, if and only if Hfifn ~5'||(A;) — )■ as n — )■ oo for any k in N. Therefore condition 
(13. ip clearly implies continuity of L. 

To establish the converse implication, note that by |Rul Theorem 1.37], a neighborhood 
basis at for C^(R) is given by 

f/fc,, = {he C^{R) I < e}, (A; G N, e > 0). 

Thus, if L is continuous, there exists for any A; in N an m in N and a positive 6, such that 
L{Um,s) ^ Uk^i- For any non-zero function g in C^(M), we have that G Um,s, and 

therefore 

|||'5|l^ilM^^II(fc) < 1' i-^- ll^^ll(fe) < III^IIm, 
which establishes (13. ip . ■ 



"(OO 



3.3 Remark. Appealing to Proposition 12. 4[ we may define a mapping S: C^(M) — C'^ 
by setting, for g in C^(]R), Sg = f, where / is the unique solution to (12. lip . By uniqueness 
S is automatically a linear mapping. We define next the linear mapping T : (M) — )■ C'l 
by the formula: 

Tg = {Sg)"\ {g G 
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3.4 Proposition. The linear mappings S,T: C^(M) — )■ C^(M) introduced in Reniark \3.3\ 
are continuous when (M) is viewed as a Frechet space as described in Definition \3.1[ 

Proof. Since differentiation is clearly a continuous mapping from C^(M) into itself, it follows 
immediately that T is continuous, if S is. 

To prove that S is continuous, it suffices to show that the graph of 5" is closed in C^(R) x 
C^(]R) equipped with the product topology (cf. |Rul Theorem 2.15]). So let (gn) be a 
sequence of functions in C^(R), such that (gn^Sgn) — )■ (g, f) in C^(]R) x C^(]R) for some 
functions f,g in C^(]R). In particular then, 

gn — g, Sgn /, and (Sgn)' f' uniformly on M as n — )■ oo. 

It follows that for any t in M, 

g{t) = lim g^{t) = lim [ gn{s)VI^ ds + - A){Sgn)'{t) + ?>t{Sgn){t)) 
= 1 r g{s)VI^^ds + {e - A) fit) + 3tfit). 

^TT J_2 



Therefore, by uniqueness of solutions to (12.111) . Sg = /, and the graph of S is closed. ■ 

3.5 Theorem. Consider the spectral density hn for G\JE{n, ^) and the hnear operator T : C^(M.) — )■ 
C^(]R) introduced in Remark \3.3[ Then for any function g in C^(M) we have that 

[ g{t)hn{t) dt = ^ [ g{t)^A-t^ dt + ^ / Tg{t) ■ hn{t) dt. 
Jr J-2 ^ Jr 

Proof. Consider a fixed function g from C^(M), and then put / = Sg, where S is the linear 
operator introduced in Remark [3.31 Recall that 

g{t) = i- j'^ g{s)V^^ds + it' - 4)/'(t) + 3tf{t), (t G M). (3.2) 

By Proposition 12.21 and partial integration it follows that 

= / fit) [n''h':it) + (4 - t')KXt) + th^it)] dt 
Jr 

= -n-^ [ f"\t)h^{t)dt- I A[/(t)(4-t2)]/i„(t)dt+ / tf{t)h^{t)dt 
Jr Jr cir 

[ - n~^r{t) - (4 - e)f{t) + 3t/(t)] h^{t) dt. 



so that 



f [{t'' -A)f\t) + ?,tf{t)]K{t)dt = \ [ f"'(t)h^{t)dt = ^ [ Tg{t)-K{t)dt. 
Jr Jr ^ Jr 
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Using f l3.2p and the fact that /i„(t) dt is a probabihty measure, we conclude that 

g{t)-K{t)dt = ^ f g{t)VA^(it+ j [{e - A)f{t) + ?,tf{t)\K{t)(it 



= g^t)VA^dt + \ [ Tg{t)-hn{t)dt, 

J-2 Jr 

which is the desired expression. 



As an easy corollary of Proposition 13.51 we may now derive (in the GUE case) Ercolani's 
and McLaughlin's asymptotic expansion \EM\ Theorem 1.4]. 



3.6 Corollary. Let T: C^{M.) — )■ C^(R) be the linear mapping introduced in Remark\K 
Then for any k in N and g in (M) , we have: 

g{t)K{t) dt = ^J2^ f [T'gmV^^dt + 4^ / [T'^g]{t) ■ h^{t) dt. 

k-1 



i=0 "^-2 



Proof. The first equality in the corollary follows immediately by successive applications of 
Theorem 13.51 To show the second one, it remains to establish that for any A; in N 



sup / \[T^g]{t) \ ■ hn{t)dt < oo. 
neN Jr 



But this follows immediately from the fact that T^'g is bounded, and the fact that hn{t) dt is 
a probability measure for each n. m 

4 Asymptotic expansion for the Cauchy transform 

For a GUE(n, ^) random matrix we consider now the Cauchy transform given by 
G„(A)=E{tr4(Al„-X„)-^]}= / -J—hnit)dt, (AgC\R). 

Setting 

g^{t) = g{\,t) = j^^, {teR, XeC\R), 

we have by the usual theorem on differentiation under the integral sign (for analytical func- 
tions) that Gn is analytical on C \ M with derivatives 

:^"(^) = I (['^t^lM t) dt = £ {§-Mt)) K{t) dt (4.1) 



dA^' 

for any A; in N and A in C \ 
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4.1 Lemma. The Cauchy transform Gn of a GUE(r;,, ^) random matrix X„ satisfies tie 
following differential equation: 

n-'^GM + (4 - A2):^G„(A) + AG„(A) = 2, (4.2) 
dA'^ dA 

for all XinC\ R. 



Proof. From Proposition 12.21 and partial integration we obtain for fixed A in C \ R that 

= / gx{t) [n-^h':{t) + (4 - t^)h'^{t) + tK{t)] dt 
Jr 

= [[- n-'gm - (4 - e)g',{t) + 2tg^{t)\ K{t) dt 

JR 



(4.3) 



Note here that 



and that 



2n / . N 4-^2 4- a2 2A 
(4-*VW=(^=(^ + ^-l. 

Inserting this into (14. 3 p and using (14. ip and the fact that /i„ is a probabihty density, we find 
that ^ 

= ^''^^^(A) + (4 - A2) ^^^(A) + AG„(A) - 2, 

as desired. ■ 

For each fixed A in C \ M, we apply next Corollary 13.61 to the function gx and obtain for 
any k in Nq the expansion: 

= /. * = *(A) + ^ + ^ + . . . + ^ + 0(„---), (4.4) 

where rij{\) = ^ J^^ [T^gx] (t) a/4 — dt for all j. To determine these coefficients we shall 
insert the expansion (14. 4p into the differential equation (14. 2 p in order to obtain differential 
equations for the ?7j's. To make this rigorous, we need first to establish analyticity of the rjj^s 
as functions of A. 

4.2 Lemma. (i) For any k in No the mapping A t— )■ T'^gx is analytical as a mapping from 
C\R into the Frechet space C^{M.), and 

d^ / 



T'gx = T' [g^giK •)) for any j in N. 



dAJ 

(a) For any k, n in N, consider the mappings 7]k, Rk,n '■ C \ M — t- C given by 

77fc(A) = [T''gx]{s)V^^ds, (A G C \ M) (4.5) 

Rk,nW = [ [T'+'gx]{s)h^{s) ds, (A G C \ R). (4.6) 
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These mappings are analytical on C\M. with derivatives: 



dX3 



%(A) = j'^[T''{di_g^X,.))]{s)V^^ds, (AgC\M, jGN) 



-Rfc,n(A) 



[T'^' {&9{\ ■))] {s)h^{s) d., (A G C \ M, J G N). 



Proof, (i) By standard methods it follows that for any A in C \ M and Z, j in No, 



lim ( 



sup 



0. 



(4.7) 



given by 



When j = 0, formula (14. 7p shows that the mapping F: C \ M — )■ 

F(A)=^?(A,-), (AgC\M), 

is analytical on C \ M with derivative ^F(A) = ^g{X,-) (cf. [Ruj Definition 3.30]). Using 
then (14. 7p and induction on j, it follows that moreover 



d^' 
dA^ 



F(A) 



9{\ 



(A G C \ M) 



for all j in N. For each in N the mapping : C^(M) — )■ C^(]R) is linear and continuous 
(cf. Proposition 13. 4p . and it follows therefore immediately that the composed mapping T'^ o 
F: C \ M — )■ (M) is again analytical on C \ M with derivatives 



d^' 



,,,T'=oF(A) = T'=(— ,(A,.)) for all, in M. 



dAJ 

This establishes (i). 

(ii) As an immediate consequence of (i), for each fixed s in M the mapping A i— t- 
\T^g{\^ ■y\{s) is analytical with derivatives 



d^' 
dA^ 



giX,-) (s), (jGN). 



Note here that by Lemma 13.21 



<c(m |-<,(A 



(m(k,0)) 



for suitable constants C(/c, 0) in (0, oo) and m(A;, 0) in N. Hence, for any closed ball B inside 
C \ R and any j in M we have that 



sup 

Ae-B 



< C{k, 0) sup 

AG-B 



9{\ 



(m{fc,0)) 



< OO. 



It follows now by application of the usual theorem on differentiation under the integral sign, 
that for any finite Borel-measure ^ on M, the mapping A t-j- j^^T^gi^X, ■)]('5) yu(ds) is analytical 
on C \ M with derivatives 



dA^ 



[T'g{\■m^^{ds) 



d^' 



-[T'^g{\■)]{s)^^{ds)= / [T\§jg{\.))]{s) ^^{ds). 



In particular this implies (ii). 
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4.3 Lemma. Let Gn denote the Cauchy-transform of hn{x)dx, and consider for each A in 
C \ R and k in No the asymptotic expansion: 



n 



(4. 



given by Corollary \3.b\ Then the coefficients rjj{X) are analytical as functions of X, and they 
satisfy the following recursive system of differential equations: 



(4-A^K(A) + Ar/o(A)=2 
(A^-4W.(A)-Ar^,(A)=r^'%(A), (j G N). 



(4.9) 



Proof. For each j in No the coefficient rij{X) is given by fl4.5l) (cf. Corollary I3.6p . and hence 
Lemma 14.21 asserts that rij is analytical on C \ M. Recall also from Corollary 13.61 that the 
Q^^-2k-2^ term in fl4.8p has the form n~^'^~'^Rk n{^), where Rk^n{^) is given by fl4.6l) and is 
again an analytical function on C \ M according to Lemma 14. 2[ Inserting now (14. 8 p into the 
differential equation (14. 2p . we obtain for A in C \ M that 

2 = n-^G';^iX) + (4 - A2)G"„(A) + AG„(A) 



n 



j=0 



n '^^'r]"j\X) + n 



-2k-2 



j=0 



+ A(5^n-2^r^,(A) + n-2'=-Xn(A) 



[(4 - X'W,iX) + Xvo{\)] + E^"'' K'^i(^) + (4 - >^"HW + 



+ n-''-' [v'^X) + (4 - X')R',JX) + XR,,nW] + n-''-'RZ^{X). 

Using Lemma 14.21 we note here that for fixed k and A we have for any / in No that 
sup|^i?,,.(A)| =sup / [T>'+^{l,g{X,-))]{s)K{s)ds <\\T'^^{9_g^X,-))\\^< 

neN neN JR 

since T''+^{-^g{X, ■)) G C^{M.). Thus, letting n ^ oo in flCTj) . it follows that 

(4-A2)r7[,(A) + Ar/o(A) = 2, (AeC\M). 
and subsequently by multiplication with that 



(4.10) 



CXO, 



= J2^-'^^"W;uw + (4 - >^'HW + H(A)] 

+ n-'' [4'iX) + (4 - A^)i?l,„(A) + Ai?,,„(A)] + 
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(4.11) 



n 



-2k-2 r>i" 



Letting then n — )■ oo in fl4.1ip . we find similarly (assuming k > 1) that 

r/[,"(A) + (4-A>;(A) + Ar/i(A) = 0, 

and subsequently that 

k 

= ,(A) + (4 - X'XiX) + Ar^,(A)] 



i=2 



+ n-''+' [r^r(A) + (4 - X')R',JX) + Ai?,,„(A)] + n-"^R'l'JX). 



Continuing like this (induction), we obtain (14 .9 1) for any j in {1, 2, . . . , fc}. Since k can be 
chosen arbitrarily in N, we obtain the desired conclusion. ■ 

For any odd integer k we shall in the following use the conventions: 

, k 



(A^ - 4)1/2 ^A^l-^, ^^2_4)/c/2^ ((^2_4)i/2)'c ^^^^2) 

for any A in the region 

n := C\ [-2,2], 

and where denotes the usual main branch of the square root on C \ (— oo, 0]. We note in 
particular that 

oo, as |A| — > oo. (4-13) 
4.4 Lemma. For any r in'L \ {—3, —4} the complete solution to the differential equation: 

(A^ - 4)^/(A) - A/(A) = ^(A^ - 4)— (a e n) (4.14) 



is given by 



(r+l)(2r+l)(2r + 3) , 2(2r + l)(2r + 3)(2r + 5) ,^,.2_,.,/2 



for all X in Q, and where C is an arbitrary complex constant. 

Proof. By standard methods the complete solution to (I4.14p is given by 

/(A) = (A^ - 4) V2 1 (A^ _ 4)-3/2-|,(A2 - 4)--V2 dA, {XeQ), (4.16) 

where /(A^ — A^^^'^-^^X'^ — 4)^^^^^'^ dX denotes the class of anti-derivatives (on f2) to the 
function (A^ — 4:)^^/'^j^{X'^ — 4)^^^^/"^. Note here that by a standard calculation, 

r A2 - 4r3/2iLrA2 - A)-r-i/2 = -(2r + 1) (2r + 2) (2r + 3) A _ 4(2r + l)(2r + 3)(2r + 5)A 
^ ^ dA3^ ^ (A2-4)'-+4 (A2_4)'-+5 

Assuming that r ^ {—3, —4}, we have (since fl is connected) for k in {4, 5} that 
A(a2 _ 4)-'--^ dA = i-^(A2 - 4)-'-^+i + C, {CeC). 
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We obtain thus that 



(A2_4)-3/2 4!.(A2-4)-^-i/2dA 

_ (2r + l)(2r + 2)(2r + 3) 4(2r + l)(2r + 3)(2r + 5) 
~ 2(r + 3)(A2 -4)^+3 + 2(r + 4)(A2-4)'-+4 + 

_ (r + l)(2r + l)(2r + 3) 2(2r + l)(2r + 3)(2r + 5) 

(r + 3)(A2 -4)'^+3 + (r + 4)(A2-4)^''+4 + ' 

where C is an arbitrary constant. Inserting this expression into fl4.16p . formula (14.151) follows 
readily. ■ 

4.5 Proposition. Let Gn denote the Cauchy-transform of hn{x)dx, and consider for each 
A in C \ M and k in No the asymptotic expansion: 

G„(A) = ,„(A) + + *W + . . . + *W + 0(„---) 

given by Corollary 13. 61 Then for A in C \ M we have that 

^o(A) = ^-^(A^-4)V2, 

r7l(A) = (A2-4)-^/^ (4.18) 



and generally for j in N, rjj takes the form: 

^.(A) = EC',,.(A^-4)— V2 

r=2j 

for constants Cj^r, 2j < r < 3j — 1. Whenever j > 1, these constants satisfy the recursion 
formula: 

(2r — 3)(2r — 1) 

Cj+i,r = ^^p^ -{{r - l)Cj,r-2 + (4r - 10)C,,,_3), (2j + 2 < r < 3j + 2), (4.19) 

where for r in {2j + 2, 3j + 2} we adopt the conventions: Cj^2j-i = = Cj^^j. 

Before proceeding to the proof of Proposition 14.51 we note that for any j in No and A in 
C \ M we have by Lemma 13.21 that 

\VjW\ < II^-'S'aIIoo < C{j,0)\\gx\\mij,o), 

for suitable constants C{j, 0) in (0, oo) and m(j, 0) in N (not depending on A). In particular 
it follows that 

\r]j{ix)\^0, as X -> cx), x eR. (4.20) 
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Proof of Proposition 14.51 The function 7]q is tlie Cauciiy transform of tlie standard semi- 
circle distribution, whicli is well-known to equal the right hand side of f l4.17p (see e.g. |VDN] ). 
Now, 7]'^'{\) = -|^(A2 - 4)1/^ so by gl]) and Lemma (with r = -1), it follows that 

r],{X) = ( - 2(A2 - Ay-'/') + C{X' - 4)1/2 ^ ^^2 _ 4^-5/2 ^ ^^^2 _ 

for a suitable constant C in C Comparing (I4.20p and (14.13^ . it follows that we must have 
C = 0, which establishes fHTTSD . 

Proceeding by induction, assume that for some j in N we have established that 

r=2j 

for suitable constants C{j, r), r = 2j,2j + 1, . . . ,3j — 1. Then by (14. 9p . Lemma 14.41 and 
linearity it follows that modulo a term of the form C{X' — 4)^/^ we have that 

3i-i 3i-i 

= ^ Q^,(!±ll(*±iK?!±3).(A2 - 4) -^-5/2 _^ ^ (^^ ^ 2(2r+l)(2r+3)(2r+5) ^^^^2 _ 4^-r-7/2 
r=2j i"=2j 

3j+l 3j+2 
= 5Z Cj,-z ^''~^^^^'~f^^'~^^ (A2 |)-^-l/2 I ^ ^^^ 2(26-5)(26-3)(25-l) ^^^^2 |-)-6-l/2 

s=2i+2 s=2j+3 
_r (aZliMliMl^lQa 4\-2j-2-l/2 , ^ 2(6j-l)(6i+l)(6i+3) / ^2 4N-3j-2-l/2 

3i+i 

+ E [(^ - ^)C,,-2 + (4s - 10)fi-,.-3] (A^ - 4)— V2. 

s=2j+3 

As before (14.201) and (I4.13P imply that the neglected term Ci}? — 4) must vanish anyway. 
The resulting expression in the calculation above has the form 

3(i+i)-i 

C,+i,(A2-4)— 

where the constants Cj+i^s are immediately given by (I4.19p . whenever 2j -|- 3 < s < 3j ' + 1. 
Recalling the convention that Cj^2j-i = = C^^sj, it is easy to check that also when s = 2j + 2 
or s = 3j + 2, formula (I4.19P produces, respectively, the coefficients to (A^ — 4)^2-'~^/2 and 
(A^ — 4)~'^-'~^/2 appearing in the resulting expression above. ■ 

Using the recursion formula (I4.19p . it follows easily that 
r/2(A) = 21(A2 - 4)-9/2 + 105(A2 - 4)-"/^ 

r/3(A) = 1485(A2 - A)-^""' + 18018(A' - 4)-^^/^ ^ 5oo50(A2 - A)-^'"'. 
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We close this section by identifying the functionals 9 ^ ^ J_2[T''5f](t)-\/4 — t"^ dt as distri- 
butions (in the sense of L. Schwarts). Before stating the result, we recall that the Chebychev 
polynomials Tq, Ti, T2, . . . of the first kind are the polynomials on M determined by the rela- 
tion: 

Tfc(cos^) = cos(fc^), {9 E [0,n], k E Nq). (4.21) 
4.6 Corollary. For each j in No consider the mapping aj : C^(M) — t- C given by 

c^Aa) = ^ ljT^g]{t)VA^dt, {g E C^{R)), 

where T : (M) — > (M) is the hnear mapping introduced in Theorem \3.5[ Consider in 
addition for each k in No the mapping Ek : (M) — C given hy 

where Tq, Ti, T2, . . . are the Chebychev polynomials given by (I4.2ip . Then for any j in N, 

k=2j ^ 

where Cj^2j, C*j,2j+i5 • • • , Cj^^j^i are the constants described in Proposition \4.5i 

From Corollary 14.61 it follows in particular that aj (restricted to C^(M)) is a distribution 
supported on [—2,2] (i.e., aj{ip) = for any function ip from C^(M) such that supp(v9) fl 
[—2,2] = 0). In addition it follows from f l4.22p that aj is a distribution of order at most 
3j — 1 (cf. |Rul page 156]), and it is not hard to show that in fact the order of aj equals 
3j - 1. 

Proof of Corollary 14.61 Let j in N be given and let denote the right hand side of 
fl4.22p . Since both aj and Aj are supported on [—2,2], it suffices to show that their Stieltjes 
transforms coincide, i.e., that 

aj{gx) = A,{gx), (AgC\M), (4.23) 

where as before gx{x) = for all x in M. Since the mapping \ gx is analytical from 
C \ M into C^(]R) (cf. Lemma [4.2p . and since the linear functionals aj, Aj : C^{M.) — t- C are 
continuous, the functions A 0-j{.g\) and A t-)- Aj{gx) are analytical on C \ M. It suffices 
thus to establish fl4.23p for A in C \ M such that |A| > 2. So consider in the following a fixed 
such A. We know from Proposition 14.51 that 

3i-i 

a,{gx)=r,,{X) = Y^C,,k{\'-^)-'-K 

k=2j 

with (A^ — 4)^^^5 defined as in (I4.12p . It suffices thus to show that 

i^.(^A) = ^(A^-4)-^-^ 
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for all k in N. So let k from M be given, and recall that g\{x) = \ Yl'^oi^)^ ^'^^ ^ 
[—2, 2]. Since ^f^f < ^^"^ since the power series 



^£(£-l)---(^-r + l); 



i-r 



converges uniformly on{2;GC||2;|<j^} for any r in No, it follows that we may change the 
order of differentiation, summation and integration in the following calculation: 



ttA 



2 . ^fc °° 



'-2 

ttA 



E 

£=0 

oo 



A/ J 



dx 



^ (f - 



e=k 



{i-k)\ 



e=k 

1 

,7r 



l)...{i-k + l)i^- 



x^-1 T.(j^^^ 



x^ 



2 7^ /£\ 

x^-^'4^ilLdx 



x^ 



A" 



Using the substitution x = 2 cos 6*, ^ G (0,7r), as well as fl4.2ip and Euler's formula for cos^, 
it follows by a standard calculation that 



T, 



V4 



dx 



((p-fc)/2) ' if P e {A; + 2m I m G Nq}, 



x^ 



0, 



otherwise. 



We thus find that 



m=0 ^ ' ^ \ / rn=0 



k\ 



m=0 



k + m 



m 



-2m-2k-l 



{2k + 2m)\ 
m\{k + m)\ 



m=0 

-2k-Am 



A 



k\ 



E 

m=0 



k + m 



m 



4 \™ 
A2 



where the third equality results from a standard calculation on binomial coefficients. Recall 
now that 



•1-2) 



'-E 



m=0 



k + m 



m 



2 1 



[z G C, \z\ < 1), 



where the left hand side is formally defined as (vT" 



-2fc-l 



, with ^/'■ the usual holomorphic 



branch of the square root on C \ (— oo, 0]. We may thus conclude that 



k\ 



A^/l-^ 



-2k-l 



k\ 



(A^-4) 



where the last equality follows from (I4.12p . This completes the proof. 
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5 Asymptotic expansion for second order statistics 

In this section we shall establish asymptotic expansions, similar to Corollary 13. 6[ for co- 
variances in the form Cov{Tr„[/(X„)], Tr„[5((F„)]}, where f,g e C^(M), X„ is a GUE(n, ^) 
random matrix and Tr„ denotes the (un-normalized) trace on Mn{C). 

For complex- valued random variables Y, Z with second moments (and defined on the same 
probability space), we use the notation: 

V{r} = E{(F-E{F})2}, and Cov{Y, Z} = E{{Y -E{Y}){Z -E{Z})}. 

Note in particular that VjF} is generally not a positive number, and that Cov{y, Z} is truly 
linear in both Y and Z. 

5.1 Lemma. Let a be a positive number, and let Xn be a GUE(n, a^) random matrix. For 
any function f from C^(M) we then have that 

V{Tr„[/(X„)]} = ^ f ifix)-fiy)fij4^,^Ydxdy, (5.1) 
where the kernel ipn is given by 

Mx,y) = 2_^fj{x)^j{y) = J- , (5.2) 

j=o V ^ X y 

and the ipj 's are the Hermite functions introduced in (12. ip . 

Proof. Formula (15. ip appears in the proof of \PS\ Lemma 3] with ip given by the first equality 
in (15. 2p . The second equality in (15. 2 p is equivalent to the Christoffel-Darboux formula for 
the Hermite polynomials (see |HTF[ p. 193 formula (11)]). ■ 

5.2 Corollary. Let X„ be a GUE(n, ^) random matrix, 
(i) For any function f from (M) we iiave that 

V{Tr„[/(X)]}= / (M^M]'p„(3;,^)dxdi/, 

where the kernel pn is given by 

Pn{x,y) = J[M^/^x)ifn~l{y^y)-^n-l{^/^x)^n{^/^y)Y■ (5.3) 
(a) For any functions f and g from C^{M.) we have that 

Cov{Tr„[/(X0],Tr4^(X.)]}= / (M^M) ^) dx d^/. 

'J^2\x-yJ\x-yJ 
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Proof, (i) This follows from Lemma [5. II by a straightforward calculation, setting = - in 

(ii) Using (i) on the functions f + g and f — g we find that 
Cov{Tr„[/(X„)],Tr„[(7(X„)]} 

= ^(v{Tr„[/(X„) + - V{Tr„[/(X„) -^7(X„)]} 

1 f {if + 9)i^) - (/ + 9)iy)T - {if - 9){x) -if - g){y)f 

■Pnix,y) dx dy 



4 Jr2 {x - y)2 

1 f Af{x)g{x)+Af{y)g{y)-Af{x)g{y)-Af{y)g{x 



4 Jr2 {x - 

fix) - fiy)\(gix) - giy) 



X — y / \ X — y 



pnix,y) dx dy, 



as desired. 



In order to establish the desired asymptotic expansion of Cov{Tr„[/(X„)], Tr„[(7(X„)]}, we 
are lead by Corollary I5.2( ii) to study the asymptotic behavior, as n — )■ oo, of the probability 
measures dxdy. As a first step, it is instructive to note that p„(x, ?/) dxdy converges 

weakly, as n — )■ oo, to the probability measure pix,y) dxdy, where 

1 4 — xy 

Pix,y) = , 2 /. 9 /. l(-2,2)(a^)l(-2,2)(l/)- (5.4) 

We shall give a short proof of this fact in Proposition 15.111 below. It implies in particular 
that if (X„) is a sequence of random matrices, such that X„ ~ GUE(n, for all n, then 

limCov{Tr„[/(X„)],Tr„[^(X„)]}= [ (^^^^^^^) (^^^^^^)p(x, ?/) dx d?/, 
n^oo 'j^^\x-yJ\x-yJ 

for all /,^eCr(M). 

The key point in the approach given below is to express the density p„ in terms of the 
spectral density hn of GUE(n, ^) (see Proposition 15.41 below). 

5.3 Lemma. Consider the functions — t- M and M — i- M given by 

1 2 

Cnix,y) = - [<^„ (x) v9„- 1 (y )-</?„-! (x)(^„(y)] , {{x,y) G M^), 

and 

n-l 

Pnix) = J2^jixf, (xGM), 

i=o 

witi ipo,fi,'P2,--- the Hermite functions given in (12. ip . We then have 

Cnix, y) = fnix)fniy) - gnix)gniy) - knix)kniy), ((x, y) G M^), 
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where 



fn{x) = liMxf + ipn~l{xf) = ^{I3nix) - x(3'^ix)), (5.5) 

9n{x) = '^{M^f - ^n-i{xY) = (5.6) 

kn{x) = V?n_l(x)(^„(x) = -^P'nix) (5.7) 

for all X inW. 

Proof. Note first that with fn,gn and /c„ defined by the leftmost equalities in f l5.5p -( 13T|l we 
have that 

fn{x) + gn{x) = (Pn{xf and fn{x) - gn{x) = ipn-l{xf, 

for all X in R. Therefore, 

1 2 

Cnix,y) = -[(^„(x)v3„_l(y) - ^n-l{x)^n{y)] 

= ^[ifn{x) + gn{x)){fniy) - ^„(i/)) + (/„(x) - ^n(x) ) (/„ (?/) + ^„(?/)) - 2A;„ (x) A;„ (y)] 
= fn{,x)fn{y) - gn{,x)gn{y) - knix)kn{y), 

for any {x,y) in M^. It remains thus to establish the three rightmost equalities in (I5.5p -( [5l7|) . 
For this we use the well-known formulas (cf. e.g. |HTlt formulas (2.3)-(2.6)]): 



^nix) = y^ifn^lix) - ^J^^ifn+lix) , {5.i 



X(pnix) = ^J^ipn+l{x) + ^ipn-l{x), (5.9) 



1 n— 1 

d 



(j^'^k{xf^ = -V2nipn{x)ipn-i{x), (5.10) 



dx 

fc=0 



which hold for all n in Nq, when we adopt the convention: = 0. 

The second equality in (15. 7p is an immediate consequence of fl5.10p . Combining (15. 8p with 
(15. 9p . we note next that 



'^'n{x) = -xifn{x) + ^/2nipri-i{x) and (/)^_i(x) = xv5„_i(x) - V2n(^„(x), 
and therefore by (15.101) 

/3;'(x) = -\/2^{^'^_^{x)^n{x) + ^n-l{x)v'n{x)) 

= -V2n{xipn-l{x)ipn{x) - \PhlLpn{xf - Xif>n-\{x)^n[x) + \p2in^ n-\{xf) 
= 2n{ipn{xf - ipn-l{xf), 
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from which the second equahty in fl5.6p follows readily. Using once more fl5.8l) and fl5.9p . we 
note finally that 

n-l 
j=0 

n-l 

i=o 

n-l 

= + (j - i)vj^i{xy - u + 

j=0 

n-2 

and therefore 

which establishes the second equality in f lS.Sp . ■ 

5.4 Proposition. Let pn be the kernel given by (15. 3p and let hn be the spectral density of 
a GUE(?2, ;i) random matrix (cf. (12. 3p ). We then have 

p4x,y) = ^[h4x)h4y) - Ah'^{x)h'^{y) - {[x,y) G M^), (5.11) 

where 

hn{x) = hn{x) — xh'^{x), (x G M). 
Proof. With (n, fniQn-, kn and /3„ as in Lemma [5.31 we have that 
Pn{x,y) = ^Cn{^/^x,^y) 

= f (/n(v^^)/n(v^2/) - 9n{^/^x)g„{^y) - K{^x)K{^y)), 
and (cf. formula (12. 3p ) 

hnix) = ^(34^x). (5.13) 
Combining (I5.13P with the rightmost equalities in (I5.5p - (l5.7p . we find that 

fn{^/^x) = :^hnix), gn{^x) = ^j^K{x) , and K{^x) = -^h'^{x), 
and inserting these expressions into (I5.12p . formula (15. lip follows readily. ■ 

By C^(]R^) we denote the vector space of infinitely often different iable functions /: — )■ 
C satisfying that 

WD'lDlfW^ := sup \D'^Dlf{x,y)\< 00, 

for any k, I in Nq. Here Di and D2 denote, respectively, the partial derivatives of / with 
respect to the first and the second variable. 



(5.12) 
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^{R) given 



5.5 Lemma. Assume that f G C^{M?) and consider the mapping (pf:M.^Cf 
by 

iPf{x) = f{xr), (xGR). 
Then cpf is infinitely often differentiable from M into C^(R), and for any k in N 

^^^^(x) = [Dlf]{xr), (xGM). (5.14) 

Proof. By splitting / in its real- and imaginary parts, we may assume that / is real-valued. 
For any /c in N the function D^f is again an element of C°°(M^). Therefore, by induction, it 
suffices to prove that v'/ is differentiable with derivative given by f l5.14p (in the case k = 1). 
For this we need to establish that 



h 



(m) 



0, as /i — )■ 



for any m in N and any x in M. This amounts to showing that for fixed a; in M and / in N we 
have that 

Dy{x + h,y)-Dif{x,y) 



sup 



h 



D^,DJ{x,y) 



0, as /i -)■ 0. 



For fixed ?/ in R second order Taylor expansion for the function [D2f]{-,y) yields that 
Dlfix + h,y)- D'jix, y) = D^D^fix, y)h + \DId'J{S,, y)h\ 



for some real number ^ = ^(x, h) between x + h and x. Consequently, 



sup 
as desired 



Dy{x + h,y)-Dy{x,y) 



h 



-DlDJ{x,y) 



< '^\\DlDlf\ 



0, as /i — )■ 0, 



5.6 Corollary. Let T be the hnear mapping introduced in Remark 13. 31 and let f be a 
function from C^(R^). We then have 

(i) For any j in Nq the mapping 

ijf. x^ T^f{x,-): R^C-t 

is infinitely often differentiable with derivatives given by 

d^ 



dx'^ 



^f{x)=T^{[D'lf]{x,-)). 



(5.15) 



(a) For any j in Nq the mapping d^- : R — )• C given by 

^A^) = ^ f [TV(x,-)](t)v/4^dt, 
is a (M.)- function. Moreover, for any k in N 



(X G 



dx 



Vj{x) 



2n 



[T^{[D'[]f{x,.))]{t)V4^dt. 



(5.16) 
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Proof, (i) As in the proof of Lemma 15.51 it suffices to prove that ip-j is differentiable with 
derivative given by f l5.15p (in the case k = 1). But this follows immediately from Lemma [53| 
since ipf = o (pj^ where T: C^(R) — )■ C^(M) is a linear, continuous mapping (cf. Propo- 
sition [^1]). 

(ii) It suffices to prove that Vj is bounded and differentiable with derivative given by 
fl5.16p (in the case k = 1). To prove that Vj is differentiable with the prescribed derivative, 
it suffices, in view of (i), to establish that the mapping 

1 , 

27r J_2 

is linear and continuous. It is clearly linear, and since 

i- r g{t)V^^dt\ < \\g\U {g e C^(M)), 

J — 2 

it is also continuous. To see finally that Vj is a bounded mapping, we note that since 
: C^(M) — 7- C^(M) is continuous, there are (cf. Lemma [3.21) constants C from (0, oo) and 
m in N, such that 

||r^/(x,-)|| <C max ||D'/(x,-)|| <C max ||M/|| 
for any x in M. Therefore, 

sup |t;j(x) I < sup ||T''/(x, •) II < C _ max ||D2/||nr, ^ 



oo (YS>2\ 



1=1,. ..,m 



since / G C^{W). u 

5.7 Proposition. For any function f in C^{M.^) there exists a sequence (/3j(/))jeNo of com- 
plex numbers such that 

[ f{x,y)K{x)K{y) dxdy = J2^ + 0{n-''~') 



for any k in Nq. 

Proof. Let k in Nq be given. For fixed a; in M the function f{x, ■) belongs to C^(M^) and 
hence Corollary 13.61 asserts that 

/ f{^.y)hn{y)dy = J2^ + ^2 I [T'^'f{xr)\{t)K{t)dt, (5.17) 
where the functions t;^ : R — t- M are given by 

^A^) = ^ f [T'fi^, ■)] {t)VA^dt, {xeR, j = l,...,k). 
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As noted in the proof of Corollary 15.61 there exist constants C from (0, oo) and m in N, such 
that 

||r'=+V(a;,-)ll <C^max \\Dif\\ , (x G R). 

M J \ ^ 'Woo — Moo' ^ ' 

Hence, since /i„ is a probability density, 



Cl := sup 



i=l m 

Using now Fubini's Theorem and (I5.17P we find that 



< C max ll-Do fll < oo. 

— , , M z oo 



= Vn-2i /" t;^.(x)/i,(x)dx + 0(n-2'=-2)^ 



(5.18) 



where the 0(n~^^~^)-term is bounded by Cln^'^^^'^. According to Corollary I5.6( ii). vj G 
C^(]R) for each j in {0,1, . . . , k}, and hence another application of Corollary 13.61 yields that 

k-j 



f v,{x)K{x)dx = J2^-^ + 0{ 



for suitable complex numbers ^o(/), . . . ,^i_j{f)- Inserting these expressions into fl5.18p we 
find that 

k k—j 

-2k-2\ 



I fix, y)K{x)K{y) <\xdy = Y,[Y. 4^ + 0{n~''-^)) + 0{n 

-'^^ j=0 1=0 ^ 

= E--"(E^M/))+o( 



r=0 ^ j=0 

Thus, setting (3r{f ) = Y7j=o C-jif)^ r = 0, 1, . . . , fc, we have obtained the desired expansion. ■ 

For the proof of Theorem 15.91 below we need to extend the asymptotic expansion in 
Proposition 15.71 to a larger class of functions than C^{M?). 

5.8 Proposition. Assume that /: — )■ C is infinitely often differentiable, and polynonii- 
ally bounded in tfie sense tfiat 

\f[x,y)\<C{l + x^ + yy\ ((x,i/)gM2) 

for suitable constants C from (0, oo) and m in Nq. Tlien there exists a sequence (/3j(/))jeNo 
of complex numbers, such that 

[ f{x,y)K{x)K{y) dxdy = J2^ + 0{n-''~') 



j=0 



for any k in Nq 
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Proof. We start by choosing a function if from C^(M^), satisfying that 

• ip{x,y) G [0, 1] for all {x,y) in R^. 
. supp(/) C [-4,4] X [-4,4]. 

• if = lon [-3,3] X [-3,3]. 

We then write / = fip + f{l-ip). Since fcp e C^{R^) C C^{R^), it follows from 
Proposition 15.71 that there exists a sequence (/3j(/))jeNo of complex numbers, such that 

/ fix, yMx, y)K{x)K{y) da; = V ^ + 0{n-^'^-') 

for any in Nq. Therefore, it suffices to establish that 



f{x,y){l-^{x,y))K{x)K{y)dxdy = 0{n 

for any k in Nq. Note here that (1 — = on [—3, 3] x [—3, 3], and that for some positive 
constant C we have that 

y){l - ^{x, y))\ < C{1 + x^ + y'T < C'ix'"^ + y'^ < C'x'"'y'"^, 

for all {x,y) outside [—3,3] x [—3,3]. Therefore, 

fix, y){l - <^(x, y))K{x)K{y) dx dy < C [ x^^y^'^K{x)K{y) dx dy 

^]R2\[_3_3]x[_3,3] 

<AC' / x2V™/^n(x)/i„(i/)dyda; 



4C"(^x2'"/i„(x)dx)(^' 



y'^'^K^y) dy 



where the second estimate uses symmetry of the function {x,y) (-)■ x^^y'^^hn{x)hn{y). By 
Wigner's semi-circle law (for moments) 

lim / x'^'^K{x)dx = — [ x^'^V^^dx, 

n^oo 271 J_2 

and therefore it now suffices to show that 

POO 

J y'^"'K{y) dy = 0{n-^^-^) for any k in Nq. (5.19) 
Recall here that hn is the spectral density of a GUE(n, ^) random matrix Xn, so that 

y'-'Kiy) dy = E{tr„ [(X„)2-l(3,oo)(X„)] } = -e{ ^^(Af ))^'n(3,oo)(Af 
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where A^") < A^^^ < ■ ■ • < A^T^ are the ordered (random) eigenvalues of Xn- Since the function 
y I—). ?/^™l(3 is non- decreasing on M, it follows that 

n 

^E(^f^)""l(3.oo)(Af^) < (A(r))^-l(3,oo)(A(r)) < ||X„f-l(3.oo)(||X„||). 

Using |HSTt Proposition 6.4] it thus follows that 

y2-/i„(y)dy < E{||X„f™l(3,^)(||X„||)} < 7(2m)ne-"/2, 



for a suitable positive constant 7 (2m) (not depending on n). This clearly implies f l5.19p . and 
the proof is completed. ■ 

5.9 Theorem. Let p„ be the kernel given by (15. 3p . Then for any function f in C^(M^) 
tiere exists a sequence (/3j(/))jeNo of complex numbers such that 



for any k in Nq- 

Proof. Using Proposition 15.41 we have that 

fix,y)pnix,y)dxdy = \ [ f{x,y)hn{x)hn{y)dxdy- / f{x,y)h!^{x)h'^{y)dxdy 



fix^y)Kix)Kiy) dy, 

(5.20) 



and it suffices then to establish asymptotic expansions of the type set out in the theorem for 
each of the integrals appearing on the right hand side. 

By Fubini's Theorem and integration by parts, it follows that 

f{x,y)K{x)h'^{y)dxdy = / ——f{x,y)hn{x)hn{y)dxdy, (5.21) 

J^2 dxdy 

and since ^^f{x,y) € C^(M^), Proposition 15.71 yields an asymptotic expansion of the 
desired kind for this integral. Similarly 

fix,y)h'^{x)h'^{y)dxdy = / f{x,y)hn{x)hn{y) dxdy, (5.22) 

where g^fg^2 f{x-, y) E C^(R^), and another apphcation of Proposition 15.71 yields the desired 
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asymptotic expansion. Finally, using again Fubini's Theorem and integration by parts, 
f{x,y)Kix)hniy) dx dy 

f{x,y){hn{x) - xh'^{x)){K{y) - yh'^{y)) dx dy 

fix, y) [K{x)hniy) - xh'^{x)hn{y) - yh'^{y)hn{x) + xyh'^{x)h'^{y)\ dx dy (5.23) 
[f{x,y) + £{xf{x,y)) + §^{yf{x,y)) + ^{xyf{x,y))]hn{x)hn{y) dxdy 

[4/(x, y) + 2x£f{x, y) + 2y-§^f{x, y) + xy^f{x, y)] hn{x)K{y) dx dy. 

In the latter integral, the function inside the brackets is clearly a polynomially bounded C°°- 
function on R^, and hence Proposition 15.81 provides an asymptotic expansion of the desired 
kind. This completes the proof. ■ 

5.10 Corollary. For any functions f,g in C^(R), there exists a sequence {(3j{f,g))j(zf,i of 
complex numbers, sucli that for any k in Nq 

Cov{Tr4/(X„)],Tr4^(X0]}= / (^^^^^^)(^-^^^^^)pn{x,y)dxdy 

jR2\x~yJ\x-yJ 

k (5-24) 

Proof. The first equality in fl5.24p was established in Proposition I5.2( ii). Appealing then 
to Theorem 15. 9[ the existence of a sequence {/3j{f,g))jeno satisfying the second equality will 
follow, if we establish that the function 



x-y ■ 

f'{x), ifx = y 



lix^y. 



^fix,y)- 

belongs to C^(]R^) for any function / from C^(]R). But this follows from the formula 

Af{x,y)= [ f'{sx + {l-s)y)ds, {{x,y) eR'), 
Jo 

which together with the usual theorem on differentiation under the integral sign shows that 
A/ is a C°°-function on with derivatives given by 

A/(x,y)= / f^'+^+'\sx + il-s)y)s\l-syds, iix,y) eR') 



dx^dy^ 
for any /c, Hn Nq 
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We close this section by giving a short proof of the previously mentioned fact that the 
measures pn{x , y) dx dy converge weakly to the measure p{x,y)dxdy given by (15 ■4p . As 
indicated at the end of the introduction, this fact is well-known in the physics literature (see 
|KKP] and references therein). 

5.11 Proposition. For each n in N, let fin denote the measure on with density pn with 
respect to Lebesgue measure on M?. Then p^ is a probabihty measure on M^, and pn converges 
weakly, as n ^ oo, to the probability measure p on with density 

Pi^'V) = ^ l(-2.2)(a^)l(-2,2)(l/), 

with respect to Lebesgue measure on M^. 
Proof. We prove that 

lim [ e'"^+™Vn(2;,l/)dxdi/= / e'"^+™V(a;,2/)dxdi/ (5.25) 

for all z,w in M. Given such z and w, we apply formulas fl5.20p - (l5.23p to the case where 
f{x,y) = Q^^^+^^y^ and it follows that 

e^^^+'"'V„(x,y) dxdy 

2 

= \ I e'^''^''"'^[hn{x)K{y) - AK{^)K{y) - n-^K{x)K{y)\ d^dy 

5.26 

1 /■ 

= - / [4 + lizx + Iwjy - zwxy + Azw - n~'^z^w^] e'^'''+""^/i„(x)/in(l/) dx dy 

= \ [ [(4 + Azw - rr'^z^w'^) + 2i2;x + 2iwy - zwxy]e'''''^''"yh„{x)hn{y) dx dy. 
A Jr2 

In the case 2; = w = 0, it follows in particular that /i„ is indeed a probability measure, and 
hence, once (I5.25P has been established, so is p. By linearity the resulting expression in 
(I5.26P may be written as a linear combination of 4 integrals of tensor products (a function 
of X times a function of y). Therefore, by Fubini's Theorem and Wigner's semi-circle law, it 
follows that 

lim / e^"^+''"V„(x,i/)dxdi/ 

= 7 / + ^-^^ + + 2iw?/ - zwxy]e''''''^''"yhoo{x)h^{y) dxdy, 
A Jm? 

where hoo{x) = ^a/4 — x'^l[_2,2]{x) ■ For x in (—2,2) it is easily seen that 

—X ~ 2 

h'ooi.^) = - — „ , and h^{x) := hoo{x) - xh'^{x) = — , (5.27) 

27rv4 — X'^ 7rv4 — a;^ 
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so in particular h'^ and hoc are both £ ^-functions (with respect to Lebesgue measure). This 
enables us to perform the calculations in fl5.26p in the reversed order and with /i„ replaced 
by hoo- We may thus deduce that 

hm f e^^^+'^ypn{x,y)dxdy= [ e'^^+'-y[\h^ix)hM - h'^i^W^oiv)] dxdy. 

"-^°°JR2 J(-2,2)x(-2,2) 

(5.28) 

Finally it follows from (I5.27P and a straightforward calculation that 

lhoo{x)hoo{y) - h'^{x)h'^{y) = ^Z^y^= (5.29) 
for all x^y in (—2,2). Combining f l5.28p with f l5.29p . we have established fl5.25p . ■ 



6 Asymptotic expansion for the two-dimensional Cauchy 
transform 

In this section we study in greater detail the asymptotic expansion from Corollary I5.1UI in 
the case where /(x) = and g{x) = for A,/i in C \ R. In this setup we put 

G„(A,/i) = Cov{Tr„[(Al - X^)-^], Tr„[(/il - Xn)-']}, 

where as before X„ is a GUE(n, -) random matrix. 

Recall from Proposition 15. 1 II that lim„_j.oo Gn{)^, /i) = G{X, fi) for any A, ;U in C \ M, where 

G(A.ri= / ( '"--'"-'"-'■^>" )( '^-^'"-'^-^''> (..,)d.dy. (6.1) 

J^2\ x-y /V x-y / 

6.1 Lemma. Let Gn be the Cauchy transform of a GUE(n, ^) random matrix. Then for any 
X in C\M. we have that 

G^{\f - AG'S\? + 4G;(A) - ^G"^(A)^ = 0, 
where Gn{\) = G„(A) - AG; (A). 
Proof. For A in C \ M we put 

ir„(A) = G„(A)^ - AG'^iXr + 4G;(A) - ^G;(A)^ 
Observing that G;(A) = -AG;(A), it follows that for any A in C \ M 

KiX) = 2G„(A)G;(A) - 8G;(A)G:(A) + 4G;(A) - -^g;(A)G:(A) 
= 2G:(A) [ - AG„(A) + A^GUA) - 4GUA) + 2 - ;^G:(A)] 
= 2G;(A) [ - ;^G:(A) + (A^ - 4)G'„(A) - AG„(A) + 2] 
= 0, 
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where the last equahty follows from Lemma 14.11 We may thus conclude that Kn is constant 
on each of the two connected components of C\]R. However, for ?/ in M we have by dominated 
convergence that 



\^yG'nijy)\ = y / ^^y^^y ^nix) dx 



< / ^ ^ K{x) dx — )■ 0, 

yZ _|_ rj.^ 



as \y\ — )■ oo, and similarly Gniyy) and G"^{iy) — )■ as |y| — i- oo. It thus follows that 
Kniiy) — )■ 0, as \y\ — t- oo, y G M, and hence we must have Kn = 0, as desired. ■ 

6.2 Theorem. Let Xn be a GUE(n, ^) random matrix, and consider for X, ^ in C \ M. the 
two-dimensional Cauchy transform: 

G„(A,/i) = Cov{Tr„[(Al -X„)-i],Tr„[(/il 

(i) If X fi, we have that 

Gn{x,f^) = ^(^Tj^i^nWGM - (2g;(a) - i){2G'M - 1) + 1 - ^g:{x)g'M), 

where G'„(A) is the Cauchy transform of Xn at X, and where G„(A) = Gn{X) — AG"„(A). 

(ii) If X = n E C\R we have that 

V{Tr„[(Al - Xnr']} = Gn{X, A) = ^(A^ - 4)G'^(A)^ - ^^^'(A)^ 
with Gn{X) as in (i). 

Proof, (i) Assume that A,/i G C \ M, and that A 7^ /i. Using Corollary 15. 2( ii) we find that 
G„(A,^)=/ ( (A - - (A - - - - 



X — y / \ X — y 

1 



{X - x){^ - x){X - y){fi - y) 

1 1 



Pn{x,y) dx dx 

Pn{x,y) dx dy. 



{fi — Xy J^2 \X — X fi — xJ\X — y p — y 
Using now Proposition 15.41 and Fubini's Theorem, it follows that 

Gn(A, /i) = ^^-l^((i7„(A) - Hn{p)r - 4(/„(A) - - ^(J„(A) - Up))') , (6.2) 

where e.g. 

Hn{X) = I hn{x)dx, /„(A) = I h'^{x)dx, and J„(A) = f h'^{x)dx. 

.Its X X ./ro A X ./m X X 
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Note here that by partial integration and fl4.ip 



X — X 



hn{x) dx 



X — X 

Gn{X) + X 
GniX). 



{hn{x) — xh'^{x)) dx = Gn{X) 

1 



A 



{X-xy 



X — X 

hn{x) dx = Gn{X) — XG'„{X) 



We find similarly that 
1 



X — X 



h'^{x) dx = G'n{X), and 



A — a; 



h'nix) dx = G"(A). 



Inserting these expressions into (16. 2p . it follows that 

4(A - M)^G(A,/i) = (G„(A) - GMY - 4(G"„(A) - G'MY - UG'LW ' ^'MY 

= 'Gn{Xf - AG'S\? - ^G:{Xf] + [gM' - ^G'M' - ^G':,{^^f 

- 2G^{X)GM + 8G;(A)GU/i) + ;^G'^(A)G'^(/i) 
= -4G;(A) - AG'M - 2Gn{X)GM + 8G'^{X)G'M + ^G:{X)G'M 
where the last equality uses Lemma I6.1[ We may thus conclude that 



Gn{X,^) 



-1 



2(A-/x)2 

-1 
2(A-/iP 



GMGM + 2G"JA) + 2G'M ' "^G'^WG'M - ^G:{X)G'M 
G„(A)G„(/i) - (2G;(A) - 1){2G'M - 1) + 1 - ^G:{X)G'M), 



which completes the proof of (i). 

(ii) Proceeding as in the proof of (i) we find by application of Proposition 15.41 that 



4G„(A,A) = 4 / 

JR 



{X-xnX~y) 



-p„(x,?/) dx dy 



h„{x) 
(X-x) 



■ dx 



(A-x)^ 



(6.3) 



(A - xy 



By calculations similar to those in the proof of (i), we have here that 

which inserted into ( 16. 3p yields the formula in (ii). ■ 

6.3 Corollary. Consider the coefEcients r]j, j G Nq, in the asymptotic expansion of Gn{X) 
(cf. Proposition 14.5)) . and adopt as before the notation fjj^X) = rij{X) — Ar7^(A). 
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(i) For any distinct A, fi from C \ {0} and k in Nq we have the asymptotic expansion: 
1 



G'„(A,yU) 



where 



2(A-/x)n 



ro(A, ■ ■+^>%^+o(n-^-^ 



n 



2k 



(6.4) 



ro(A,/i) = {2i,{\) - i){2i,{^) - - 1, 

and for I in {1,2, . . . , k} 

i-i 



r,(A,/i) = 24{X){27^M - 1) + 2r/K/x)(2r/[,(A) - 1) + Aj^Vjim-^if^) 

i-i I 

j=0 j=0 

(the third term on the right hand side should be neglected, when 1 = 1). 
(a) For any X in C\M. and any k in Nq we have that 

G.(A, A) = I [To(A) + ^ + ^ + -- - + ^ + 0(n---) 
4 L n'^" 



(6.5) 



where 

and for I in {1,2,..., k} 



T,(A) = (A^ - 4) J]r/;(A)r/;V(A) -J^V-'WVi-i-.W- 

j=0 j=0 

Proof. From the asymptotic expansion of G'„(A) (cf. Proposition 14. 5p it follows that 
2G;(A) - 1 = (2r^;(A) - 1) + MP + . . . + MP + Oin-^'-'^ 



(6.6) 



n 



2k 



GniX) = ^o(A) + 



^i(A) 



+ ■■■ + 



^2k 

VkjX) 

^2k 



0{n 



-2k-2\ 



where we also use that the derivatives of the remainder terms are controlled via Lemma 14. 2[ 
Inserting the above expressions (and the corresponding expressions for 2G"„(/i) — 1, G'^{fi) 
and GnifJ')) into the formula in Theorem I6.2( i). it is straightforward to establish (i) by 
collecting the n~^'-terms for each / in {0, 1, . . . , A;}. The proof of (ii) follows similarly from 
Theorem 16.2^ 11) . ■ 
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6.4 Remark. Using that ?7o('^) = f ~ |('^^ ~ 4)^/^ (cf. Proposition I4.5P it follows from 
Corollary I6.3( i) and a straightforward calculation that for distinct A and /i from C \ M, 

G(A,/.) = jim G„(A,/.) = = ^(X^ ( (p- 4)^2 _ 4)1/2 " l)' (6-7) 

where G{X, fi) was initially given by (16. ip . If A = yU, it follows similarly from Corollary I6.3( ii) 
that 

G(A, A) = Jim G„(A, A) = ^(A^ - ^WoW' = p^^^, 

which may also be obtained by letting A tend to /x in (16. 7p ! 

Using also that = (A^ — 4)~^/^, it follows from fl6.5p and a rather tedious calculation 
that 

ri(A, fi) = _ 4)7/2 (5 V + 4A V + V - 52A/i=' + 3A=^/i=' - IQfi'' + 4A V 

- 52AV^ + 208A/i + 5A V - 52AV - ISA^ + 320 + 4A^) . 
Inserting this into (16 ■4p and letting A tend to fi we obtain that 

Ti(A) =4(21A2 + 20)(A2-4)-^ 
which is in accordance with (I6.6p . 
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